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Abstract 

A classical counterexample due to E. De Giorgi, shows that the weak maximum principle 
does not remain true for general linear elliptic differential systems. After that, there are some 
efforts to establish the weak maximum principle for special elliptic differential systems, but 
the existing works are addressing only the cases of weakly coupled systems, or almost-diagonal 
systems, or even some systems coupling in various lower order terms. In this paper, by contrast, 
we present maximum modulus estimates for weak solutions to two classes of coupled linear 
elliptic differential systems with different principal parts, under considerably mild and physically 
reasonable assumptions. The systems under consideration are strongly coupled in the second 
order terms and other lower order terms, without restrictions on the size of ratios of the different 
principal part coefficients, or on the number of equations and space variables. 
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1 Introduction 

Let m,n € N\{0}, and ft C M m be a bounded domain with an C 1 boundary L and having the 
cone property. We consider the following nonhomogeneous, isotropic elliptic differential system of 
second order: 



-div^Vy 1 ) - div(a 12 Vy 2 ) 
-div^Vy 1 ) -div(a 22 Vy 2 ) 



div(a ln Vy") + ^ C u • Vy i + D 1 ■ y = f l in ft, 

i=l 
n 

div(a 2 "Vy") + ^ C 2i ■ Vy l + D 2 ■ y = f 2 in ft, 



i=l 



(1.1) 



div(a nl Vy 1 ) - div(a" 2 Vy 



div(a" n Vy™) + s >TC ni -Vy i + D n -y = f n in ft, 



8=1 



k y 1 = g\ y 2 = g 2 , 



y n = g n 



on L, 
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and the following general nonhomogeneous elliptic differential system of second order: 

m n 

E [KdM + + • • • + + E ° u ■ v y l + 1)1 • y = f 1 in n > 

>,q=l i=l 
m n 



p,g=l i=l 

(1.2) 



- E + + • • • + w»s P )x,J + E c ™ • w + D n ■ y = r m n, 

p,q=l i=l 

I j/^s 1 , y 2 = g 2 , y n = g n onT. 

In both (1.1) and (1.2), y = (y 1 ,--- ,y n ) T is unknown, while a* J , ap ? , C ij ', D*, /* and 5* (i,j = 
1, ■ ■ ■ , n; p, q = 1, • • ■ , m) are suitable given functions (See the next section for the assumptions on 
these functions). The main purpose of this paper is to study the weak maximum principle, or the 
boundedness of weak solutions, for systems (1.1) and (1.2) with suitable measurable principal part 
coefficients. 

It is well-known that the weak maximum principle is one of the basic issues in the theory of 
partial differential equations and it plays an essential role in the study of many other problems. 
For example, a central problem in the calculus of variations is the regularity of stationary points 
for functionals of the type 

J(u) = / F(x,u(x),X7u(x))dx, 
Jn 

where u(x) = (n 1 (x), u 2 (x), ■ ■ ■ ,u n (x)) T is a vector-valued function defined on S7, and F(-, •, •) is a 
suitable function defined on $7 x R n x M mn . Research in this area has been stimulated by D. Hubert's 
Problem 19, which can be reduced to the regularity of weak solutions to nonlinear elliptic equations 
or systems. This problem was successfully solved by C. Morrey ([14]) in two dimensions and the 
general case with n = 1 was finally solved by E. De Giorgi ([7]) and J. Nash ([16]), and refined 
by J. Moser ([15]). We refer to [1, 2, 9, 10, 11] and the references cited therein for more details 
in this respect. A fundamental step of the De Giorgi-Nash-Moser approach in solving the scalar 
Hilbert's 19th problem (i.e., n = 1) is to establish the weak maximum principle for single linear 
elliptic equations. 

In many physical and geometrical applications, u may be a vector function, and therefore, 
the corresponding Euler-Lagrange equation is a system. Naturally, one expects to extend the De 
Giorgi-Nash-Moser approach to the case of systems. However, in 1968, E. De Giorgi ([8]) gave a 
surprising counterexample of an unbounded solution to a second order linear elliptic system with 
bounded coefficients. This means that the weak maximum principle fails for general second order 
linear elliptic systems, and therefore, the De Giorgi-Nash-Moser estimates are no longer valid for 
general elliptic systems. 

In order to establish the weak maximum principle for elliptic differential systems, as a conse- 
quence of the above mentioned De Giorgi's counterexample, one has to impose some restrictions 
on the structure of the system. There exist a few works in this direction. In [11], a weak maximum 
principle was proved for a class of special elliptic systems with variable coefficients, in which the 
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principal operator in each equation takes the same form, and it is acting only on one component 
of the solution vector. In [3, 4, 12], some weak maximum principles were discussed in the frame of 
Campanato's space for linear or quasilinear elliptic systems under some additional conditions, say, 
2 < m < 4 in [3], the coefficients matrix being constant in [4], and a dispersion assumption on the 
eigenvalues of the principal part coefficients matrix in [12] (and hence the system is almost-diagonal 
in high space dimensions). 

In this paper, we choose the usual Sobolev space as the working space and derive weak maximum 
principles for two classes of strongly coupled elliptic systems with different principal parts, in the 
spirit of the classical framework for single equations. We emphasize that our systems are strongly 
coupled, i.e., the (second order) terms of the principal parts are coupled each other. Therefore, 
when establishing the desired a prior estimate, it is necessary to get rid of some undesired terms 
generated by different principal operators and/or different solution components appeared in the 
same equation. This goal is achieved by choosing delicately suitable weighted test functions. As 
far as we know, this is the first result on the weak maximum principle (in the classical sense) for 
strongly coupled elliptic systems. 

The rest of this paper is organized as follows. Section 2 is devoted to stating the main results 
in this work. In Section 3, we collect some preliminary results which will be useful later. Sections 
4 and 5 are addressed to the proof of the main results, i.e., the boundedness of weak solutions 
to systems (1.1) and (1.2), respectively. Finally, in Section 6, we give an example in which the 
assumptions for proving the boundedness of the weak solution to system (1.2) are satisfied. 

2 Statement of the main results 

To begin with, we introduce some assumptions. Suppose that, for i,j = l,--- ,n, 

0*3 G L°°(n) (2.1) 

and 



(2.2) 



C ij (-) e L e (n-R m ) and D i (-) € L^(Q;W l ) for some 6 > m, 
f = (f\ ■■ , P) T e H-\n- IT), g = {g\ ■ ■ ■ , g n ) T € H\Q; R n ), 
and for i,j = 1, ■ ■ ■ , n and p, q = 1, ■ ■ ■ , m, 

a%EL°°(n), u'! q a',',. (2.3) 

Moreover, we assume that, for some positive constant p, 

n m 

i,j=i P =i (2.4) 



and 



»,J=lp,3=l (2-5) 
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Conditions (2.4) and (2.5) mean that both systems (1.1) and (1.2) are elliptic (see [5, Section 1 
of Chapter 8]). Clearly, system (1.1) is a special case of system (1.2). The weak solution to system 
(1.2) is understood in the following sense: 

Definition 2.1 We call y = (y 1 , • • • , y") T € i^ 1 (f2; IR") to be a weak solution to system (1.2) if for 
any<p = {<p 1 ,--- ,<p n ) T € H%(Sl;R n ), 



n m 



„ n n 

E E / 4{(x)yi pV i q dx + / ^[E^w-w+^w-^ 

i,j=lp,q=l JSl Jn i=l j=l 



dx 



— (/>^}tf-l(f};R™),tfl(n;]R»); 

and y % — g % € Hq(Q), i = 1, • • ■ , n. 

Similar to the proof of [5, Theorem 2.3 in Chapter 1]), it is easy to show the following well- 
posedness result for system (1.2). 

Lemma 2.1 Let conditions (2.2), (2.3) and (2.5) be fulfilled. Then, there exists a constant vq = 
vo(n,m,9) > such that system (1.2) admits a unique weak solution y € 1 ($7; IR") whenever the 
following inequality 



26 

e-m I ,2 



E (D l (x) ■ M ) // > [ E |C^I^(f2;R-) 

V (x, //) = (x, fj, 1 , [i 2 , ■ ■ ■ ,fi n ) 6fix 



i=l 



(2.6) 



ZioWs for v > uq . Moreover, 

\y\m(n-M") < c(n,m,n,p, \a l J q \ La o^, \C tj \L<>(n-SL™)> \ D% \ L % m-R«)) (l/l-ff _1 («;R n ) + Mff 1 ^™))- 

The proof of Lemma 2.1 is standard and therefore we omit the details. 

Remark 2.1 Clearly, if C l;, (-) = for i,j = 1, ■ ■ ■ ,n, i/ien condition (2.6) is satisfied whenever 
the function matrix (D 1 (x), • • • ,D n (x)) is semi-positive definite. 



Next, we put 



A 



( a 11 a 21 
a 12 a 22 



y a ln a 2n 



a" 2 



/ a 22 a 32 
a 23 a 33 



y a 2n a 3 " 



a™ 2 \ 

a n3 



/ 



Also, denote by B lJ (i,j = 1, • • • ,n) the cofactor of ^4 with respect to a 13 and by detyl the 
determinant of matrix A. It is easy to see that B 11 = B. Moreover, under condition (2.4), it is 
easy to show that det B / 0. 

The first main result in this paper is the following boundedness of weak solutions to system 
(1.1). 
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Theorem 2.1 Suppose that conditions (2.1), (2.2) and (2.4) are fulfilled, inequality (2.6) holds for 
v>vq (Recall Lemma 2.1 for vq), f £ L2(fi;IR n ) and 



det B ij 
detB 



tW 1 * 00 ^), i,j = l,---,n. 



(2.7) 



Then the weak solution y € ^(fi; M n ) to system (1.1) satisfies 

esssup \y\ < C I m, n, 9, tt, p, \a tJ \ L ^ (n) , \C lJ \ L e (n . Rm) , \D l \ L <> , 



det B l 3 



det B 



\g\m m , l/l L S (n . Rn) > esssup 2/| 
The proof of Theorem 2.1 will be given in Section 4. 

Remark 2.2 We conjecture that assumption (2.7) in Theorem 2.1 is a technical condition, and 
therefore it is not really necessary. However, we do not know how to drop this assumption at this 
moment. 

Since almost all of the natural materials are isotropic, Theorem 2.1 suffices for most of physical 
applications. Nevertheless, from the mathematical point of view, it would be quite interesting to 
extend Theorem 2.1 to more general anisotropic systems such as (1.2), in which the scalar functions 
a%3 (hj = 1, 2, ■ ■ ■ , n) appeared in system (1.1) are replaced by the W nxm matrix-valued functions 
(%>?) i<p q < m - Note however that, by the above mentioned De Giorgi's counterexample ([8]), this 
seems to be highly nontrivial in the general case. In the rest of this section, we shall extend Theorem 
2.1 to system (1.2) under some technical assumptions. 

In order to treat system (1.2), we put 

\ 



Mpg = 



We assume that 



/ 


a 11 

u pq 


a 21 ■ 






a 12 

u pq 


a 22 ■ 

u pq 


■ a n2 

u pq 



V a 



In 
pq 



In 
l pq 



L 



pq 



det M, 



pq: 



(p,q = !,■■■ ,m). 



Lpq 7^ 0, 



(2.8) 



Also, we denote by v pq = 1, ■ ■ ■ , n; p, q = 1, ■ ■ ■ , m) the cofactor of M pq with respect to a pq . 
Further, let us introduce the following assumption: 

(H) There exist functions f pq , h ij G W 1 ' 00 ^) (i, j = 1, • • ■ , n; p, q = 1, • • • , m) such that 
1) h 11 = 1, h 1 ^ = h? 1 , and the following matrix is uniformly positive definite: 



V 



( 1 h 12 
h 21 h 22 



\ h nl h n2 

i.e., V > pilnxn for some positive number pi; 
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h ln \ 

h 2n 

h nn J 



f k 

2) The function E l i := y~] h^v 1 ^ is independent of p and q, and E 1 ^ G W 1 ' 00 ^) ft 

L pq , =1 



or any 



i,j = !,■■■ ,»v 

3) T/ie following matrix is uniformly positive definite: 

( F\\ F12 ■ ■ ■ F\m \ 

F 2 l F22 • • • Fi m 



\ F m i F m 2 



Fmm / 



i.e., F > P2l rnxm for some positive number pi, where F pq 
1, • • • ,m; 

4) The following matrix is uniformly positive definite: 

( F h 12 F ■ ■ ■ h ln F \ 

h 21 F h 22p ... h 2np 

M 



E 

i=i 



a pq El1 f° r an V P> Q 



i.e., M > P3I; 



nmxnm 



\ h nl F h n2 F ■ ■ ■ h nn F J 
for some positive number p^. 



Now, we can state our another main result as the following boundedness result for weak solutions 
to system (1.2). 

Theorem 2.2 Suppose that conditions (2.3), (2.5) and (2.8) are fulfilled, C^(-) G L°°(^;M m ), 
£)*(•) g L°°(r2;M n ) and inequality (2.6) holds for v > u ((Recall Lemma 2.1 for Vq)), f G 
L°°(r2;]R n ), and assumption (H) holds. Then the weak solution y G H l (Q,;M. n ) to system (1.2) 
satisfies 

esssup \y\ < C (m, n, Q, p, p 1 , p 2 , p 3 , |a^| L oo( n ), |C lJ | i0O ( n . R m), |L> l | L oo (Q . Mn) , 



P'liyi.^n) ' l^lw 1 . 00 ^)) bln 1 ^")) l/U°°(r2;R") 5 esssup |y| J . 

The proof of Theorem 2.2 will be given in Section 5. Also, in Section 6, we shall give an 
illustrative example, in which all of the assumptions in Theorem 2.2 are satisfied. 

Remark 2.3 It is well-known that one of the classical topics in partial differential equations is the 
strong maximum principle for elliptic differential equations, which has many applications ([6, 17, 
18, 19] and so on). However, the existing results on strong maximum principle are mainly focusing 
on single elliptic equations, although one can find some works on weakly coupled elliptic systems 
([1, 13, 20]) and the references therein. It would be quite interesting to establish a strong maximum 
principle for system (1.1) or even for system (1.2), but this remains to be done and it seems to be 
far from easy. 
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3 Some preliminaries 

In this section, we collect some known preliminary results which will be useful later. 
The first one is the following interpolation result. 

Lemma 3.1 ([11, Theorem 2.1 in Chapter 2]) For any u G Wj'*(fi), t>l and r > 1, it holds that 

\u\ LP * {n) < ^IVul^njIuliTj,)' 

where a 



^7 — — pr) 1 , t* = and (3 is a constant depending only on m, t, p* , r and 

a. Moreover, if t < m, p* can be any number between r and t* ; if t > m, p* can be any number 
larger than r. 

For any Lebesgue measurable function u defined on Q, we put = {x G £l;u(x) > k} and 
denote by {A^l the Lebesgue measure of set A^. The next lemma is quite useful in deriving the 
supremum of function u. 

Lemma 3.2 ([11, Theorem 5.1 in Chapter 2]) Suppose that u G W l > ma {VL) n L qo (n) for some 

tuq G [l,m] and some qo > 1. If for any fixed k > esssupu, function u satisfies the following 

r 

inequality: 



/ |V«pdx < 7 f (u-k) 

J A k \_J A k 



+ 1 k°\A k \ 1 -^ +s \ 



where 7, Iq, a and eo are positive constants satisfying lo < and mo < cr < eoqo + mo, then 

esssup-u < C*(fi, m , q , 7, '0, o - , £o, esssupu, |n| LTO rn))- 
n r 

Moreover, when a = mo, |«|i,9o(n) appeared in C* can be replaced by |«|x,i(n)- 

The last lemma is a result on comparison of the determinants between a matrix and its sym- 
metrizing matrix. 

E + E T 

Lemma 3.3 ([21, Theorem 3.7.1]) For a real matrix E, if H(E) = is positive definite, 

then 

det H(E) < det E. 

4 Proof of Theorem 2.1 

The goal of this section is to prove our first main result, i.e., Theorem 2.1. 

Proof of Theorem 2.1. First, for the weak solution y = (y 1 ,?/ 2 ,-- - ,y n ) T to system (1.1), any 

fixed k > esssup \y\ 2 and r > 0, put 
r 

4> r (x) = mm{(\y(x)\ 2 - k) + ,r}, 
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where s+ = maxjs, 0} (for any s G M). We choose (f = (if 1 , if 2 , • • • , </? n ) T as a test function, where 

ip* = (y 1 ^ +y 2 T 2i -\ Yy n T ni )(\) r and T ij = 1, 2, ■ ■ ■ ,n) are suitable functions to be specified 

later. Obviously, <p € i?o(fi;R n ). By Definition 2.1, it follows that 

n „ 

V / {a ll Vy l + a i2 Vy 2 + ■■■ + a in Vy n ) ■ V [(y 1 ^ 1 + y 2 T 2% + ■■■ + y n T ni )<p r ] dx 

n „ n 

+ E / ( E ° ij ■ v ^ + ° l ■ y) (y lTlt + y 2T2i + • • • + y" Tm )^^ 
i=i ^ j=i 

n „ 

This implies that, for any E\ > 0, 

/ { {a^Wl 2 + a i2 T 2i \Vy 
i=i ^ n 



+ •• - + a m T ra |Vy n | 2 ) 0, 



+ ^ (a ij T li Vy j ■ Vy l <j> r + a ij T li y l Vy j ■ V0 r ) 

J,j€{l,2,-,n},/^j 

+1 [ a <i T ii V ( y i)2 + a i2 T 2i v(y 2 )2 + . . . + a ^T ni V(y n ) 2 ] • V^ r }dx 

n . n 

E / [/* (^ Tli + y 2 ? 21 + ■■■ + y nTni )<t>r - E °V vj/' • vT^> r 
i=i i,j=i 

n 

- ( £ • Vy^' + D* ■ y) (y 1 T li + y 2 T 2t + ■■■ + y n T m )4 



(4.1) 



<C I \f\\y\<t> r dx + £l [ \Vy\ 2 c/> r dx + C 
In Jn Jn 



i=i 



|y| 2 (/) r (ix. 



(4.2) 



Here and hereafter C denotes a generic constant (which may be different from one place to another), 
depending only on n, m, p, |a y |^oo(Q) and |T* J |wi,oo(n) (z, j = 1,2,--- , n). Also, it is clear that 

[a il T H V(y 1 ) 2 + a i2 T 2i V(y 2 ) 2 + ■■■ + a in T ni V (y 11 ) 2 ] ■ V<p r 

= o^T^Vlyl 2 • V0 r + [(a i2 T 2i - o^T^JV^ 2 ) 2 + • • • + (a in T ni - a a T u )V (y n ) 2 ] ■ V<p r . 

Next, we choose suitable functions T % i = 1,2, ••• , n) such that T* J 6 W 1,00 (fi) and some 
undesired terms in (4.1) vanish. For this purpose, we consider the following linear system: 

n 

JV'T" = 0, V j, I = 1,2, • • • ,n with j / I, 

i=l 
n 

J^(a i2 T 2i - a jl T") = 0, 



(4.3) 



J2(a in T ni - a il T u ) = 0. 



v i=l 



By Lemma 3.3 and (2.4), it follows that detyl > p n and deti? > p n 1 . One can easily check that 
the following is a solution to system (4.3): 



det B li 

r n = 1 rii = ( _i)i+<^LL_ (i = 2> . 

det B 

ri' = (-i)'+'y ffi ¥^ (j = 2,3,- 



,n) 



(4.4) 



Z=i 



Also, it is not difficult to check that 



,n, i = 1,2,--- ,n). 



det A 
det B' 



detB 3% 

From this fact and noting (4.4), we see that T JJ = (— 1) J+JI — - — — (i,j = 1,2, ■■■ ,n). Moreover, 

d.et -D 

there exists a constant p* > 0, depending only on n,p and |a u |l<»(q) = 2, • • • ,n), such that 



n 

Inrpnl \ r 



^ a il T lI = . . . = ^ a «n T n« > 



i=l 1=1 

Combining (4.2), (4.3) and (4.5) with (4.1), we arrive at 



det B 



>P*- 



(4.5) 



/ (\Vy\ 2 cp r + |V<M 2 ) dx < C [ \f\\y\(f> r dx + / [ V |C^| 2 + V |^| + 1 | |y| 2 <Mx 



• (4-6) 



In the sequel, we estimate the terms in the right side of (4.6) and prove that the left side of 
this inequality is uniformly bounded with respect to r > 0. First of all, by Holder's inequality and 
Lemma 3.1, we see that 



/ \f\Wrdx< I \f\(\y\ 2 -k^^dx + k^ [ 

Jn Jn Jn 

(\y\ 2 -k)hr r 44, +Ci|/| 



< / 49 

L^(Q;R») 

< l/l 49 

L^+6(Q;R») 

< l/l 49 

L3+6(f};R") 



2 



(4.7) 



M 2 - ^)<Ar 



+ 1 



'L^=2(f2) 

here and hereafter C\ denotes a constant depending only on k, 0, n, m and $7. Put L 

n n 

y~] [C* J ; |i9(Q. R m) + ^ 1-^*1 r £^ m ^ + 1- Lemma 3.1, we find that 



i=l 



'L2(n ; R") 



pin n 



|y| 2 <?V< 



v i,j=l i=l 

< dL / (|y| 2 r )<^cte 
<CiL|(|y| 2 -A;)0 r | e 



<CiL|(|y| 2 -fc)</»J « +CiL|^| 9 



(4.8) 



+ c , i- L lylii(n ; i;n ) - 
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On the other hand, put = \/(|y| 2 — k)(j) r . It follows that 



/ \Vu*\ 2 dx = / \X7u*\ 2 dx = 
Jn J\y\ 2 >k Jk 



< 2 



\y\ 2 >k 



\y\ 2 >k 



'\y\ 2 >k 



2^ r yVy + (\y\ 2 -k)V<f> r 



2V(|y| 2 - k)^ 



dx 



dx + 



2^ 



2 Jk+r>\y\ 2 >k 



k)V(j)r 



dx 



<2 / (\y\ 2 -k)- 1 <p r \y\ 2 \Vy\ 2 dx + 2 [ 

J\v\ 2 >k Jk 



'\y\ 2 >k 

Noting 4> r < \y\ 2 — k, we see that 



k+r>\y\ 2 >k 



(\y\ 2 -k)^ l \y\ 2 \Vy\ 2 dx. 



(\y\ 2 -ky 1 <i )r \y\ 2 \Vy\ 2 dx= / <t) r \Vy\ 2 dx + k (|y| 2 - k)' 1 <t> r \V y\ 2 dx 

\y\ 2 >k J\y\ 2 >k J\y\ 2 >k 

< / 4> r \Vy\ 2 dx + k i \Vy\ 2 dx< / r |Vy| 2 dx + / \Vy\ 2 dx. 
Jn J\y\ 2 >k Jn Jn 

Noting that <f> r = \y\ 2 — k whenever k + r > |y| 2 , it is clear that 

/ (\y\ 2 -k)^ l \y\ 2 \Vy\ 2 dx = / \y\ 2 \Vy\ 2 dx 

Jk+r>\y\ 2 >k Jk+r>\y\ 2 >k 



-L 

Ik 



k+r>\y\ 2 >k 



-k)\Vy\ 2 dx + k [ |Vy| 2 dx 

Jk+r>\y\ 2 >k 



\Vy\ 2 dx + k / \Vy\ 2 dx< / (f) r \V y\ 2 dx + k \ \Vy\ 2 dx. 

k+r>\y\ 2 >k Jk+r>\y\ 2 >k Jn Jn 

Therefore, by (4.9)-(4.11), we conclude that 

f |Vu*| 2 cix<4 f |Vy| 2 r dx + Ci f \Vy\ 2 dx. 
Jn Jn Jn 

By (4.12) and Lemma 3.1, for any < £2 < 1, we end up with 

\{\y\ 2 -k)(j) r \ T _e_ f ^= (^J ut^dx^j - S2 J Nu^dx + CiEz 1 (^J \u*\dx^ 



'L#=2(fi) 



f(\y\ 2 -k)hl 
Jn 



dx 



<4e 2 / \Vy\ 2 cl) r dx + C 1 e2\y\ 2 H i {n . R n ) + C 1 £ 2 
< 4e 2 J \Vy\ 2 (j) r dx + C^ly^i {n . Rn) + C , ieJ 1 |y|i 2(n . R „ ) . 
Therefore, substituting (4.13) into (4.7) and (4.8) respectively, we see that 

4e 2 ^ |Vy[ 2 r dx + Cie2|y|Hi(f7;R") + Ci £ 2 1 bll2(Q;iR«) " 
+Cl|/ y (f W y| ^(^) 



»l^< \f\ L & {a , 
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and 



~ / n n \ 

/ £i^'i 2 +£i^i + 1 

JQ \ij=l i=l / 



\y\ 2 (f) r dx 



< ±C x Le 2 J \Vy\ 2 <f> r dx + C\L(1 + Ci^)^ 2 ^.^ + C 2 L£ 2 l \y\ 4 L2{n . Rn) . 
Combining the above inequalities with (4.6) and taking e 2 sufficiently small such that 



4 l/l 4^ e 2 + 4CiLe 2 C < -, 

L5+6(f7;R") / 2 

where C and Ci are the constants appeared in (4.6) and (4.8) respectively, we arrive at 

/ {\Vy\ 2 4>r + \V4>r\ 2 )dx<C 2 , 

JQ 

here and hereafter C 2 is a constant depending on C, C\, L, \ f\ e and \y\ii 1 (Q-R n )> independent 

of r. Since 4> r G Hq(Q), by the definition of (f> r , letting r — >■ +00 in the above inequality, for any 

fixed k > esssup |y| 2 , we obtain that 
r 

[ \Vy\ 2 (\y\ 2 -k) + dx+ f[(\y\ 2 -k) + ] 2 dx+ [ \V(\y\ 2 )\ 2 dx < C 2 . (4.14) 
JQ JQ J\y\ 2 >k 

Finally, we construct a sequence of inequalities with respect to A k , where A k = {x G O; |y(x)| 2 > 
k}. Again, by (4.6), we get that 

f (|Vy| 2 r + |V0 r | 2 ] dx 

JQ 

I \f\\y\(\y\ 2 -k) + dx + I j^i^C 13 ? + \D l \ + l)\y\ 2 {\y\ 2 -k) + dx 

JQ JQ^l^f^ ' 



< C 



Letting r — > +00 in the above inequality, for any e 3 > 0, by the Holder inequality and Lemma 3.1, 
we see that 



/ |Vy| 2 (|y| 2 -fc)dx+ f \V\y\ 2 \ 2 dx 

J A k J A k 

„ re ra 

<C [|/| + ^(^|^| 2 + |^| + l)j(|y| 4 + l)dx 



<C(\f\ e +L) ( \\y\ 2 -k 



'L?(Q;WL n ) 



2 20 + k 2 \A k \ 1 - 2 e ) +C(\f\ +L)\A k \ 1 - 2 o ( 415 ) 



< H' + L + l/l , _ . , ) (e 3 |V| y | 2 || 2(Afc) + C(e 3 )\\y\ 2 - k\\ 2(Ak) ) 



+C ( + e J fc 2 |A fc 



2| A. |l-4 
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Denote v = \y\ and take £3 to be sufficiently small, then by (4.14) and (4.15), one derives that 



f \Vv\ 2 dx<C 3 [ \v-k\ 2 dx + C 3 k 2 \A k a - 
J A k J A k 



(4.16) 



where C3 denotes a constant only depending only on C, L and \ f\ L ^ nRn y 
By Lemma 3.2, we take 

2 2 

to = 2, l = 2, cr = 2, e = 7 = CV 

to y 



Then it follows that 

esssup|y| < C to, n, 9, tt, p, |a lj | L ->(n), L, 



det 



detB 



M^ 1 '°°(S1) 



(4.17) 



I^L^R")' eSS ^p|2/|) 

Since y is the weak solution, by Lemma 2.1, we have that 

ML 2 (n;R") ^ C ( m > n ' fi > P> l aU U°°(fi)> L )(l/l L | (n . B;n) + k!//l(n ; R"))- 

This, combined with (4.17), yields the desired conclusion in Theorem 2.1. □ 

5 Proof of Theorem 2.2 

Now, let us prove our second main result, i.e., Theorem 2.2. 

Proof of Theorem 2.2. The main idea is the same as that in the proof of Theorem 2.1. First, 
for any weak solution y = (y 1 , y 2 , ■ ■ ■ , y n ) T to system (1.2), We choose (p = ((p 1 ,^ 2 , • • • , ip n ) T € 

n 

#rj(S7;R n ) as a test function, where (p l = '^^E tl y l ( r , and E' l i = 1,2, ■■■ , n) are given by 

1=1 

assumption (H), while Cr is a suitable function to be specified later. By Definition 2.1, it follows 
that 

n m n „ n „ n „ 

E E E/ apX(^VCrM*+ E / ^-V^U + E / D'-yEUj/trdx 

i,j=lp,q=ll=l Jn i,j,l=l Jfl i,l=l Jfl 

= E / fE a y l ( r dx. 
This implies that 



n m n 



»J=lp,5=l '=1 



E E E / *^ p l4,Cr + CI^l4y , (Cr)x,+aa(^),,wi p y , Cr 



n n ,. 

+ E / ^'-V^VCrdx+E / D'-yE^Crdx 

n ,. 

V / fE u y l ( r dx. 



12 



Therefore, 



n m n 



»,j=lp,5=l '=1 

< c 4 



/ 



l/IMCr + (i+E 1^1 1 \Vy\\v\Cr + Y.\ Dl \\y\ 2 tr 

*J=1 



(5.1) 



dx, 



here and hereafter C4 denotes a constant depending only on n, m, p, {apgl^oo^ and \E 13 '\w l >°°(Si) 
= !,■■■ ,n;p,q = 1, • • ■ ,m). 

Next, we estimate the terms in the left side of (5.1). For this purpose, by (2.8), condition 
2) in assumption (H) and the Cramer rule, we see that for any p,q = l, -- , m, functions E 13 

n 

(i, j = 1, • • • , n) (given by assumption (H)) satisfy E o,p q E lj = f pq h l 3 . In particular, by h 11 = 1, 

n 

we see that f pq = E^ U - Therefore, 



1=1 



1=1 



Y j a l ; q E l3 = h i3 Y J ^ q E 11 . 



(5.2) 



1=1 



1=1 



Hence, 



TL Tfl TL n lib iv n / iv \ 

EEE ai 3 q E il yi p y l (C r ) Xq dx= £ £ / E<^ UW^kA 

i,j=lp,9=l Z=l ■ /n p,g=li,/=l- /n \i=l / 

= E /JE"^ 1 ) 

P,9=l ,/S2 \i=l / 
m „ _. / n 

-E E«: 



l£^V') 2 + E WW 

i =1 j,ie{l,2,-,n}, j</ 



(Cr)^dx (5.3) 



./n 2 . 



E(^ 7 MV 
3,1=1 



(Cr)x q dx. 



On the other hand, by condition 4) in assumption (H) and noting (5.2), it is easy to see that 
/ 



M 



E«n^ ••• E«n^ 



1=1 



1=1 



\ 

Ell rpln » 
a lm E 

1=1 



E«n^ E 



In T-i/1 



E 

1=1 



An jpln 



E<^' 



in 



a ml E 
\ i=i 



E 



In rpll 
u mm J - J 



EAn rpln 
a ml E 

1=1 



E a mm E ' 



In 



1=1 
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Therefore, 



n m n ,. ,. 

E E E / a%E il yi p y l Xg C r dx > p 3 \Vy\ 2 ( r dx. 
.,j=lp,q=l l=i Jn Jn 



(5.4) 



Combining (5.3) and (5.4) with (5.1), we have 



1/lMCr + [ i + E l^'l ) \ v y\\y\(r + E l^lfol'Cr + bi 2 |vcr 



<c 5 



»J=1 



1=1 



(5.5) 



where ip = h? l yiy l and C5 depends only on C4, ps and ; |wi,°°(n) j = 1, 2, • • ■ , n). 
i,i=i 

For s, r > and A; > sup -0 s , denote 

r 

A k = {x G Q, I V s (a:) > ^} , A£ = {a; € O | k < i/j s (x) < k + r} . 

Moreover, we choose ( r = min{r, — k) + }. Then, by (5.5), and using condition 3) in assumption 
(H), we conclude that 



I \Vy\ 2 ( r dx + [ i) s ~ l \Vi)\ 2 dx 

JA k JAV 



>A k 
<C 6 



l/IMCr + ( 1 + E \° l '\ 2 + E 1^1 I \y\^ + bPlVCr 
»J=1 



i=l 



(5.6) 



dx, 



where Cq denotes a constant depending only on s, C5 and pi- Moreover, using condition 1) in 
assumption (H), it follows that 

ip > pi\y\ 2 - 

Now, let us estimate the right side of (5.6). First, by Holder's inequality, for any £4 > 0, we 

have that f zti f 

/ \f\\y\Crdx < C 7 / \y\( r dx<e A ( r >dx + e?C 7 \y\ s+1 dx, (5.7) 

J A k J A k J A k J A k 

n 

here and hereafter C7 denotes aconstant depending only on Cq, Q, pi, \f\L°°(n-R n ), E \L oc (Q;R m ) 

n 

and ^ \ d1, \l°°(SI$l»)- Next, 
i=i 

f 11+ E l^f + El^l ] M 2 Crd*<C 7 / |y| 2 C r dx 
■ /j4 * \ ij=i i=i / ^ 

<e 4 / CT^dx + e^C-j f \y\ 2s+2 dx. 
J A h JA k 



(5.8) 
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Further, 



f \y\ 2 \V( r \dx <C 7 [ \y\ 2 r~ l \VMdx<eJ' ^{V^dx + £j X C 7 f |y|V _1 ^ 
J A k Jai Jai Jai 

£4/ ^{V^dxdx + e^C 7 [ \y\ 2s+2 dx. 

J AT J AT 



< 

On the other hand, by Poincare's inequality, 

2 



f Cr'dx <C 7 [ 
Jn Jn 



3+1 

V(Cr 2s ) 



dx<C 7 I Cr s V s ~ 2 |VVf ^ < C 7 I ^ s - l \V^\ 2 dx. (5.10) 



Therefore, by (5.6)-(5.10), taking £4 sufficiently small, we get that 



/ \Vy\ 2 Crdx+ [ \V(^)\ 2 dx<C 7 I (\y\ 2s+2 + \y\ s+1 )dx. 
JA h JAI JA k 



>A k JA r k JA k 

Letting r — > +00 in the above inequality, we have that 

[ \Vy\ 2 {r ~k)dx+ f \V{^)\ 2 dx<C 7 f (\y\ 2s+2 + \y\ s+1 )dx. (5.11) 
J A k J A k J A k 

2 

Notice that for any given constant s < (if m < 2, s can be any positive number), the right 

m — 2 

side of (5.11) is finite. 

~ s-f-1 s+1 ~ 

Denote k = k~ and = {x G | ip~ > k} = A k . Then by (5.11), it follows that 
[ \V(^)\ 2 dx < C 7 [ 4> s+1 dx + C 7 [ ip^dx 

<C 7 [ -~k) 2 dx + C 7 k 2 \A~ k \ + C 7 ! (^ - ~k)dx + C 7 ~k\A~A (5.12) 

<C 7 [ - kfdx + C 7 k 2 \A~ k \. 

By Lemma 3.2, we take 



s+i ~ 2 

u = Y> 2 f rno = a = /o = 2, k = k, 7 = C7, £0 = — • 

to 

Then, using also Lemma 2.1, it follows that 

esssup \y\ < C (to, n, Q, p, pi, p 2 , P3, \a% q \L°°{Q), \C tj \ L °°(n;Rm), |-D l |z,°°(n ; R»), 



I^IwloCO)' l^lwi.-(f2)>lflifi(n ; ]R«), |/|L<»(n;R»). esssup|y|). 
This completes the proof of Theorem 2.2. 
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6 An example 



In this section we give an example, in which the coefficients a v 3 q = 1, 2, 
of system (1.2) satisfy all of the assumptions in Theorem 2.2. 

For any given functions b %3 € VF 1,00 (f2) and g pq £ L°°(ti) = 1,2, 
such that g pq > and the following matrix is uniformly positive definite: 

/ 911 912 ••• gim \ 



,n;p,q = 1,2, ■ 
,n;p,q = 1,2, • 



,m) 
,m) 



G :-- 



921 922 



92r, 



we take 



b ij g Pq , 



h ij 



V 9ml 9rrt2 ' ' ' 9mm 

I 1 if « = j 



'P9 



Then it is easy to check the following assertions: 

i) Condition 1) in assumption (H) holds, since V = I n xn', 

ii) By the definition of a pq and f pq , and b %] € W l, °°(£l) (i, j = 1, 2, • • • , n;p, q = 1, 2, • • ■ , m), if 
(2.8) holds, condition 2) in assumption (H) is satisfied; 

„ n 

iii) If b 11 / 0, b a = 0, i = 2,3,- ■■ ,n, and £ n := He > in 17, then by (2.5), 

conditions 3) and 4) in assumption (H) hold true. Notice that by the definition of E 11 , 
E 11 > if and only if 

/ b 22 b 23 ■■■ b 2n \ 

^32 ^33 ... yin 



det 



> 0. 



(6.1) 



y yn2 ^n3 . . . ynn j 

Moreover, under condition (6.1), the hypothesis (2.8) also holds; 
iv) Condition (2.5) is equivalent to that the following matrix is uniformly positive definite: 

\ 



/ b u G 

-b 12 G 
2 



K 



lb l2 G 



-b ln G 
2 



b 22 G 



1. 



1 



b"G -(b S2 + b 26 )G 



_( fo 2n + 6 n2 )G 



(ft3n + 6 n3) G 



1 



1 



b Ln G _ { J b 2n + h n^ G 



b nn G 
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Notice that if for some constant > 0, 

6" > p* and b ij < np* (i, j = 1, 2, • • ■ ,n; i ^ j), (6.2) 

then the matrix K is uniformly positive definite. 

By the above assertions i)-iv), suppose that the coefficients Op q = 1,2, ••• , n; p, g = 
1, 2, • • • , m) of system (1.2) satisfy that 

a pq = b J g pq , 

where 6« G W 1 ' 00 ^), c/pq € L°°(fi), 6 11 ^ 0, ft* 1 = (i = 2, • • • ,n), 5pg > 0, and G is uniformly 
positive definite, and (6.1)-(6.2) are satisfied. Then, by Theorem 2.2, we conclude the boundedness 
of weak solutions to the corresponding system (1.2). 
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